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Lattice Approach to Diquark Condensation in Dense Matter 

Susan Morrison (UKQCD collaboration) a * 

a Department of Physics, University of Wales Swansea, 
Singleton Park, Swansea SA2 8PP, U.K. 

Results are presented of a Monte Carlo simulation of a three dimensional Gross-Neveu model with SU(2) ® 
SU(2) chiral symmetry at non-zero baryon chemical potential fi, corresponding to non-zero baryon density. The 
phenomenon of quark pair condensation is sought via measurement of a two point function. For /i sufficiently 
large there is a sharp transition between a phase where the chiral symmetry is broken by a condensate {qq) and 
one where a scalar diquark condensate \ {qq) \ 7^ 0. Global U(l) baryon number symmetry may remain unbroken, 
however, due to the absence of long range order in the phase of (qq). 



Recent work ]T],|| suggests that new forms of 
ordering may exist at high baryon density. In ad- 
dition to the standard (qq) condensate associated 
with spontaneous breaking of chiral symmetry it 
has been proposed that the existence of a Fermi 
surface combined with any attractive qq interac- 
tion results in the formation of a {qq) conden- 
sate of quark Cooper pairs. Remarkably, this di- 
quark condensate is associated with baryon num- 
ber violation. Formation of a diquark condensate, 
which is non- invariant under the vector like U(l) 
global symmetry of baryon number is generally 
forbidden by the Vafa-Witten (VW) theorem |J. 
Notable exceptions to the VW theorem include 
theories with complex measures and models with 
Yukawa couplings to scalar degrees of freedom, 
such as the coupling to a in ([j]) below. 

For /j, 7^ the fermion determinant of QCD 
is complex and lattice simulations fail dramati- 
cally however the lattice approach can be suc- 
cessfully applied at ^ to toy field theories. 
I will present work done in collaboration with 
S.J. Hands which is the first lattice study of di- 
quark condensation at fi 7^ implemented in the 
Gross-Neveu (GN) model with SU(2)®SU(2) chi- 
ral symmetry and formulated in d = 2 + 1 space- 
time dimensions. With the exception of confine- 
ment this model incorporates all of the essential 
physics. 

The lattice model we have simulated has the 



following Euclidean path integral: 



Z = J DaDTfdet(M f M[a, vr] 



>e s 2 (1) 



where M is the staggered fermion matrix; a and 
the triplet 7? are real auxiliary fields defined on 
the dual sites x of a three dimensional Euclidean 
lattice. The lattice model is formulated such that 
the integration measure is real and positive. The 
full lattice action (involving isodoublet staggered 
lattice fermion fields x an d C) an d its symmetries 
are fully discussed in j| . The corresponding d = 
3 + 1 model is outlined in Q . 

In the limit m — > 0, our model is invariant 
under a global symmetry akin to the SU(2) L ® 
SU(2) JJ of the continuum NJL model. For T = 
ji = 0, a mean field treatment shows that this 
SU(2) <E> SU(2) symmetry is spontaneously bro- 



ken to SU(2) i 



by the generation of a conden- 
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sate £0 = ( <T ) = -p ixx) an d the fermion acquires 
physical mass mf = So -The symmetry is broken 
for coupling g 2 > g 2 ~ 1.0. The numerical re- 
sults presented here are for a value g 2 = 2.0, cor- 
responding to a theory deep in the broken phase 
for ji — 0.0, and with So = 0.706(1) in units of 
inverse lattice spacing, ie. rather far from the 
continuum limit. 

For T, fi 7^ a mean field solution is also known 
|q], in which S(/x, T) is expressed in terms of So- 
At T = the basic feature is that S remains 
constant as \i is increased up to a critical value 
fjt c = So, whereupon S falls sharply to zero (in the 



2 



chiral limit), signalling a first order chiral symme- 
try restoring transition. Monte Carlo simulations 
H support this picture, even when massless Gold- 
stone excitations (with the quantum numbers of 
the 7r held) are present Jt]]. 

We performed simulations of the model (|l|) on 
L 2 x Lt lattices with bare mass to = 0.01 at a 
coupling 1/g 2 — 0.5. The simulation method is 
a standard hybrid Monte Carlo algorithm. We 
monitored the chiral condensate (xx)> an d the 
baryon density n using stochastic estimators. 

We have constructed two local diquark op- 
erators which are anti-symmetric in spatial, 
flavour and isospin indices. Firstly the spec- 
tral scalar: X p ( x ) T 2° q X q ( x ) which is even under 
parity and secondly the spectral pseudo-scalar: 
(—) xl+X2+X3 X p ( x ) T 2 q X 9 ( x ) which is odd under 
parity. 

Measurement of diquark condensates on the 
lattice demands innovative techniques and for this 
initial study we chose an indirect method: the two 
point function (qq(0)qq(x)) was measured as the 
expectation of the product of two fermion prop- 
agators as in QCD baryon spectroscopy, and the 
condensate extracted by assuming clustering at 
large spacetime separation, ie. 



{qq{0)qq{ x )) = (qq(Q)qq( x ))c + (qq)(qq), 



(2) 



the latter term on the right being proportional to 
\(qq)\ 2 . A non-zero condensate signal therefore 
shows up as a plateau in the timeslice correlator 
G(t) at large time separation, and may be ex- 
tracted by a fit of the form 

G{t) = ^2(qq(0)qq(t,x 1 ,x 2 )) 

= Aexp{~M + t) + Bexp(-M-(L t -t)) 
+C 2 (L s )\(qq)\ 2 7 (3) 

with M± the masses of forward and backward 
propagating diquark states respectively. The 
value of the constant C(L S ) is not determined a 
priori. Naively one expects C = L s , and hence 
the plateau height to be an extensive quantity. 
The trend in the diquark condensates as we pass 
from broken to symmetric phase is clear from 
Fig. |. There is a critical value, fi c ~ 0.65, 
at which the chiral (qq) condensate falls sharply 
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Figure 1. Overview of observables (xx)j n an d 
diquark condensate signals C|(<7g}| as functions 
of n, obtained on a 16 2 x 24 lattice, showing the 
onset of diquark condensation at /i c ~ 0.65. 

and the number density (n) begins to rise from 
zero. The spectral scalar C\(qq)\ condensate rises 
slowly from zero for 0.4 < /i < /i c , jumps up- 
wards discontinuously close to /i c , and continues 
to rise steadily. For \i > 1.0 as the number den- 
sity approached saturation the scalar qq propaga- 
tors became distorted. The spectral pseudoscalar 
C\(qq} \ signal was consistent with zero for [i < [i c 
and no fits to the propagators were possible but 
for fi > /j, c this condensate is non-zero and in- 
creases with /i. It is considerably smaller in 
magnitude than the spectral scalar. This pseu- 
doscalar condensate is parity violating, and there- 
fore it would be remarkable if the signal in the chi- 
rally symmetric phase were to remain non-zero in 
the large volume limit. 

The trend in the M± data |j was clear with 
small masses in the symmetric phase and large 
masses in the broken phase. The forward prop- 
agating M + states were more difficult to extract 
from the fits than the backward propagating M_ 
states. For /i < fx c , M_ decreased linearly with 
fi reflecting the trend observed in the physical 
fermion mass to/ in previous simulations [Q . M— 
was approximately constant for /i > /i c . The zero 
density fermion mass to/ = So — 0.7, defines the 
ratio of physical to cutoff scales and the observed 
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Figure 2. Spectral scalar diquark correlator at 
/it = 0.8 for both free and interacting quarks on 
L 2 S x 24 lattices. 

qq states had masses around 0.3 which is compar- 
atively light. 

Fits to the spectral scalar correlator in the sym- 
metric phase (/j, = 0.8) were very stable as L t was 
increased from 16 to 40. The smaller parity vio- 
lating signal, on the other hand, decreased with 
increasing L t . The non-zero spectral scalar con- 
densate in the broken phase is at least in part 
due to finite lattice size. To determine the be- 
haviour with spatial volume we next performed 
a series of runs on L 2 x 24 lattices at /j, = 0.8, 
with L s varying from 8 to 24. To our surprise 
we found very little change as L s increased, the 
fitted value of C\(qq)\ more or less saturating for 
L s > 16. Results for L s — 8, 16, 24 are shown 
in Fig. |2|, together with the fit for L s = 24. For 
comparison we also plot the equivalent correlators 
for free fermions at fj, — 0.8; it is striking that for 
this case the trend as L s increases is in the oppo- 
site direction. We conclude that the constant C 
is roughly independent of L s 

The fact that the two-point correlation function 
exhibits clustering at large temporal separation is 
suggestive of diquark condensation at fi ^ 0. In 
order to quantify the measurement, however, it 
will be necessary to provide a numerical estimate 
for the constant C in (|J). 

The spontaneous breaking of the U(l) sym- 



metry usually implies the existence of a Gold- 
stone mode associated with long wavelength fluc- 
tuations in the phase of the {qq) condensate. 
However, in the absence of an explicit diquark 
source (analogous to a bare mass for the case of 
(qq)), one would expect large finite volume ef- 
fects, generically proportional to L~( d ~ 2 \ where 
d! is the dimension of the effective field theory 
describing the fluctuations of the order parame- 
ter ||. We therefore speculate that the effective 
field theory describing the spatial correlations of 
the diquark correlator has d! = d—1 = 2, and that 
there is no long range order in the spatial direc- 
tions A similar distinction between temporal 
and spatial correlations has been observed in in- 
stanton liquid models Q . This would account for 
the volume-independence of C. Strictly speaking, 
therefore, we have not observed a true condensa- 
tion. This interpretation of our results can be 
tested by simulations of the equivalent 3 + 1 di- 
mensional model H , where we would expect long 
range order. An extension of this work will in- 
volve the introduction of a diquark source term 
to study the one point function. This will serve 
to calibrate our two point function measurement 
and to extend our study to other diquark opera- 
tors and the T^0 regime. 
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